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RESUMEN 


Los hongos son los principales degradadores de la madera en los ecosistemas bos- 
cosos, contribuyendo significativamente al ciclo global del carbono. Los enfoques 
metagenómicos basados en un amplicón específico (metabarcoding) constituyen una 
herramienta poderosa para su prospección y estudio. El objetivo principal de este 
estudio fue caracterizar a través de secuenciación de segunda generación (NGS), 
comunidades de hongos degradadores en la albura de dos especies de Nothofagus 
de los bosques del norte de la Patagonia, para evaluar patrones de diversidad en 
sitios, estaciones, hospedadores, compartimentos del árbol y condiciones de salud, 
como una contribución a la autoecología de las especies de este género. Nuestro 
estudio comprendió tres pasos metodológicos principales: (1) muestreo de madera 
de árboles vivos sanos y enfermos de las especies forestales N. dombeyi y N. pumi- 
lio; (ii) extracción de ADN, amplificación y secuenciación de la región ITS1 en la 
plataforma MiSeq lllumina, (iii) procesamiento de lecturas y extracción de datos de 
los órdenes Polyporales e Hymenochaetales y (iv) análisis de datos e interpretación. 
Se obtuvieron un total de 35 unidades taxonómicas (variantes de secuencia de am- 
plicón -ASV-), las cuales fueron asignadas a 23 géneros de hongos putativos en 15 
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familias. Postia pelliculosa fue la especie detectada con mayor frecuencia en el estudio. 
El hospedador fue el factor más fuerte entre las variables estudiadas en cuanto a su 
efecto sobre la estructura y la composición de la comunidad fúngica analizada. Para 
N. dombeyi, que se distribuye en una amplia gama de condiciones climáticas, el sitio 
fue el modelador más fuerte de sus comunidades, mientras que para N. pumilio se 
observó una mayor susceptibilidad a los cambios de temperatura y estacionalidad, 
que son, ciertamente, factores relevantes para la conservación de los bosques en el 
actual escenario de cambio climático. Este es el primer estudio que utiliza NGS 
como una estrategia rápida y a gran escala para descubrir la diversidad de hongos 
que degradan la madera en los bosques templados de Patagonia. 


Palabras clave — Decaimiento; Hymenochaetales; Illumina; Nothofagus; Polyporales. 


ABSTRACT 


Fungi are the main degraders of wood in forest ecosystems, contributing signifi- 
cantly to the global carbon cycle. Metagenomic approaches based on a specific am- 
plicon (metabarcoding) constitute a powerful tool for their prospecting and study. 
The main objective of this study was to characterize, through second-generation 
sequencing NGS, communities of degrading fungi in the sapwood of two species of 
Nothofagus from the forests of northern Patagonia, to evaluate diversity patterns in 
sites, seasons, hosts, plant compartments and health conditions, as a contribution 
to the autoecology of the species of this genus. Our study comprised three main 
methodological steps: (1) sampling of wood from healthy and diseased living trees 
of the forest species N. dombeyit and N. pumilio; (ii) DNA extraction, amplification 
and sequencing of the ITS1 region on the MiSeq Illumina platform, (iii) read pro- 
cessing and data extraction of the orders Polyporales and Hymenochaetales and (iv) 
data analysis and interpretation. A total of 35 molecular operational taxonomic units 
(amplicon sequence variant -ASV-) were obtained, to 23 putative fungal genera in 
15 families. Postia pelliculosa was the most frequently detected species in the study. 
The host was the strongest factor among the variables studied in terms of its effect 
on the structure and composition of the fungal community analyzed. For N. dombeyi, 
which is distributed in a wide range of climatic conditions, the site was the strongest 
shaper of its communities, while for N. pumilio a greater susceptibility to changes 
in temperature and seasonality was discovered, which are, in fact, relevant factors 
for the conservation of forests in the current scenario of climate change. This is the 
first study to use NGS as a rapid, large-scale strategy to elucidate the diversity of 
wood-degrading fungi in Patagonian temperate forests. 


Keywords — Decay; Hymenochaetales; Illumina; Nothofagus; Polyporales. 
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INTRODUCCIÓN 


Los bosques andino patagónicos son una de las últimas reservas mundiales de bos- 
ques templados con poca alteración antrópica y valiosa diversidad (Arroyo et al., 
1996). Entre los 20 géneros endémicos de árboles presentes, el 90% del área total 
está ocupada por especies de Nothofagus Blume (Donoso Zegers, 1993). Nothofagus 
dombeyi es la especie dominante en los bosques templado-húmedos de la Patagonia 
Norte y Nothofagus pumilio domina los rodales de altura en la Patagonia Norte y Cen- 
tral (Daniels y Veblen, 2004; Donoso Zegers, 1993). Ambas especies de Nothofagus 
presentan una situación distintiva de muerte en pie progresiva que exhibe un patrón 
espacial agrupado y del cual se han recuperado patógenos fúngicos potenciales, pero 
sin conducir a ninguna etiología patógena fúngica primaria (de Errasti et al., 2015; 
Molina et al., 2020; Pildain et al., 2009). Hay antecedentes, acumulados durante las 
últimas décadas de investigación, sobre la diversidad de hongos asociados a estos 
árboles en ciertos grupos, como políporos, ectomicorrizas, hongos ofiostomatoides 
y otros hongos de descomposición (Barroetaveña et al., 2019; de Errasti et al., 2016; 
Rajchenberg, 2006). Se sabe que esta micobiota juega un papel clave en la aptitud y 
el funcionamiento de los árboles a través de dinámicas complejas (Baldrian, 2016) 
y expresadas a lo largo de un continuo de mutualismo, comensalismo y parasitismo 
que puede transcurrir, incluso, durante la vida del mismo organismo fúngico (Ro- 
binson et al., 2004). En este sistema, los hongos degradadores de madera (HDM) son 
patógenos relevantes y reconocidos para los árboles muertos y en crecimiento. 

Los HDM son capaces de descomponer componentes complejos de alto peso 
molecular en la pared celular de la planta en pequeñas moléculas para su asimilación 
por los hongos involucrados, así como por otros microorganismos (Blanchette, 1991). 
Juegan un rol esencial en el reciclado de nutrientes, ya que la descomposición de la 
madera es un proceso ambiental sumamente crítico y representa el principal recurso 
de carbono en los ecosistemas y ejerce una gran influencia en el desarrollo de los 
suelos (Blanchette, 1991). La mayoría de los hongos son saprótrofos que viven en 
madera muerta. Sin embargo, algunos son patógenos, que pueden infectar árboles 
vivos y causar pudrición o descomposición en los tallos y raíces de los árboles. Los 
HDM durante mucho tiempo han sido categorizados como hongos de pudrición 
marrón, blanca y blanda en base a las preferencias de sustratos y el patrón de de- 
gradación. Esta separación es arbitraria y actualmente se conoce que los tipos de 
pudrición son más un continuo (Riley et al., 2014), pero la categorización es útil 
para la clasificación general de estos agendes degradadores y para su discusión. Los 
hongos de pudrición blanca exhiben excelentes habilidades de descomposición y 
son los responsables de la degradación de la lignina, celulosa y hemicelulosa (Tuor 
et al., 1995), contribuyen al 90% de los hongos que causan la descomposición de la 
madera y son muy diversos, particularmente en áreas tropicales (Dai, 2005). Los 
hongos de la podredumbre parda se caracterizan por ser rápidos degradadores de 
celulosa y hemicelulosa, y acceden a los polisacáridos de la planta a través de una 
potencial modificación o degradación de la lignina (Hatakka, 2005). Estos hongos 
son los mayores invasores de la biomasa forestal, representando el 10% de los hongos 
que causan la descomposición de la madera, distribuyéndose principalmente en zonas 
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templadas y frías. Varios estudios han reportado que los hongos de la pudrición parda 
han evolucionado de un hongo saprofítico de pudrición blanca (Arantes y Goodell, 
2014). Tanto los hongos de pudrición blanca como los hongos de pudrición parda son 
Basidiomicetes. Los hongos de pudrición blanda son principalmente Ascomicetes 
que en presencia de excesiva humedad descomponen la celulosa y las hemicelulosas 
y no se considera un problema en los árboles vivos en pie. 

En las últimas décadas, las nuevas tecnologías de secuenciación de alto rendi- 
miento (HTS, por sus siglas en inglés) se volvieron frecuentes y accesibles, lo que le- 
vó al progreso en las investigaciones de micobiomas vegetales, especialmente a través 
de enfoques de metabarcoding. De todos modos, los esfuerzos en la investigación de 
micobiomas y los antecedentes disponibles en dicho tópico se distribuyen de manera 
desigual entre los distintos compartimentos vegetales, las regiones geográficas y los 
biomas (Harrison y Griffin, 2020). En lo que respecta a los biomas forestales, se han 
realizado varios estudios en hoja, mientras que la micobiota degradadora de madera 
cuenta hasta la fecha con muy pocos estudios en ecosistemas de bosques templados 
del hemisferio sur. En este contexto el objetivo de este trabajo fue caracterizar la 
micobiota degradadora de la madera pertenecientes a los órdenes Hymenochaetales 
y Polyporales de Basidiomycota en los bosques de Nothofagus describiendo sus pa- 
trones de diversidad y detectando las variables de estado sanitario, estación del año 
y de hospedador que afectan sus estructuras. 


MATERIALES Y MÉTODOS 
Área de estudio y muestreo 


El estudio se realizó en el Parque Nacional Los Alerces (PNLA), Argentina (42° 58’ 
27.075” S 71° 38 37.725” O), de noviembre de 2017 a mayo de 2018. El muestreo se 
repitió estacionalmente al principio y al final de la temporada de crecimiento (de- 
nominada primavera y otoño, respectivamente) para las especies de Nothofagus. Se 
definieron como sitios de muestreo cuatro rodales de N. dombeyi y tres rodales de N. 
pumilio y, todos ellos con mortalidad progresiva en pie (Fig. 1). Se muestrearon diez 
árboles de diámetro a la altura del pecho (DAP) similares en cada sitio por mues- 
treo, por temporada: 5 sintomáticos y 5 asintomáticos según lo definido en Molina 
et al. (2020). Se muestrearon un total de 140 árboles. De cada árbol seleccionado 
se recolectaron muestras de 5 mm de diámetro y 2 cm de largo de madera (albura) 
del fuste al DAP y de las raíces utilizando un barreno de Pressler esterilizado con 
etanol al 70% (v/v) y flameado entre cada muestra. Las muestras se almacenaron 
individualmente en bolsas de plástico y se mantuvieron a -20 °C hasta su posterior 
procesamiento, un máximo de 48 horas después. 
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Fig. 1. Sitios de estudio. A) Area de estudio en el Parque Nacional Los Alerces. Sitios de muestreo 
seleccionados para los muestreos definitivos. Círculos: N. dombeyi. Triángulos: N. pumilio. MS: 
sitio río Alerce; ME: sitio lago Menéndez; AC: sitio cerro Alto el Petiso (rodal de coihue); K: sitio 
lago Kruger; AL: sitio cerro Alto el Petiso (rodal de lenga); Ri: sitio cerro El Riscoso; D: sitio cerro El 
Dedal. Mapa creado con QGIS de código abierto y gratuito. La grilla refleja coordenadas geográficas 
decimales (proyección Gauss-Krúger). B) Manchones de mortalidad agrupada en algunos sitios de 
estudio ejemplificados con flechas rojas. 


Fig. 1. Study sites. A) Study area in Los Alerces National Park. Sampling sites selected for final sam- 
pling. Circles: N. dombeyi. Triangles: N. pumilio. MS: Alerce river site; ME: Lake Menéndez site; CA: 
Alto el Petiso Mt. site (coihue stand); K: Lake Kruger site; AL: Alto el Petiso Mt. site (lenga stand); Ri: 
El Riscoso Mt. site; D: El Dedal Mt. site. Map created with QGIS. The grid shows decimal geographic 
coordinates (Gauss-Krúger projection). B) Grouped mortality patches in the study (red arrows). 


Extracción de ADN, preparación de librería 
y secuenciación de segunda generación 


Para cada muestra se molieron aproximadamente 50 mg de madera combinando la 
inmersión en nitrógeno líquido y 3 minutos de agitación en un molino mezclador 
según la metodología de Doyle y Doyle (1990), adaptada por Dumolin et al. (1995). 
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El ADN total se extrajo con el kit DNeasy Power Plant Pro (QIAGEN, Hilden, 
Alemania) siguiendo las recomendaciones del fabricante después de 10 minutos de 
incubación a 65 °C. 

La librería del espaciador interno transcrito 1 (ITS1) se preparó utilizando la 
estrategia de etiquetado dual en dos pasos TrueSeq. La amplificación de ITS1 se rea- 
lizó utilizando el par de cebadores IT'S1-F y ITS2-R (Gardes y Bruns, 1993; White et 
al., 1990) y MylaqTM Mix (Bioline, USA, Inc., Memphis) en un volumen total de 
25 uL por reacción con las siguientes condiciones de ciclcado: 94 °C durante 5 min, 
32 ciclos de 94 °C por 45 seg, 50 °C por 45 seg y 72 °C por 1 min, y una extensión 
final a 72 °C por 7 min. Los productos de PCR se purificaron utilizando ExoSap-IT 
(USB Corporation, Cleveland, OH) siguiendo las instrucciones del fabricante y se 
indizaron mediante el uso de combinaciones de barcodes específicos de la muestra 
de los pares de cebadores TruSeq 15-TS-DI-5xx e 17-TS-DI-7xx (Integrated DNA 
Technologies, Inc., Coralville, IA) con las condiciones de ciclado: 95 *C durante 3 
min, 8 ciclos de 95 °C durante 30 s, 55 °C durante 30 s y 72 °C durante 30 s, y una 
extensión final a 72 °C por 5 min. Los productos de PCR se purificaron con ExoSAP- 
IT (Thermo Fisher Scientific) siguiendo las instrucciones del fabricante y el ADN 
se cuantificó con un espectrofotómetro NanoDrop (ThermoFisher, Waltham, MA). 
Los controles negativos de PCR y extracción de ADN, así como comunidades simu- 
ladas sintéticas no biológicas SynMock (Palmer et al., 2018) (GenScript, Piscataway, 
NJ) se incluyeron como muestras y también se amplificaron, se incluyeron en los 
conjuntos finales y se secuenciaron (Palmer et al., 2018). Las muestras se agruparon 
en proporciones aproximadamente equimolares. Las bibliotecas de la región ITS1 
se secuenciaron en las instalaciones centrales de Purdue Genomics (Universidad de 
Purdue, West Lafayette, IN) con un kit de reactivos MiSeq v2 de 500 ciclos en la 
plataforma Illumina MiSeq (2 x 250 pb). 


Análisis bioinformático 


El procesamiento de los datos de salida del secuenciador se llevó a cabo con el pro- 
grama Amplicon toolkit (AMPtk) (v1.2.4; Palmer et al., 2018) que demostró ser de 
mejor rendimiento para el análisis de amplicones IT'S fúngicos en estudios previos 
con comunidades simuladas (Anslan et al., 2018; Nilsson et al., 2019) y el flujo de 
trabajo (pipeline) más adecuado para este conjunto de datos particular según Molina 
(2022). El pipeline aplicado primero eliminó las lecturas cortas (Palmer et al., 2018), 
luego, recortó las secuencias de los cebadores de las lecturas; luego, fusionó las lectu- 
ras de extremos emparejados mediante USEARCH (v9.2.64; Edgar, 2010) y realizó el 
filtrado de calidad de las lecturas ensambladas mediante el recorte de errores espera- 
dos (Edgar y Flyvbjerg, 2015). Posteriormente, constituyó las unidades taxonómicas 
(variantes de secuencia del amplicón: ASVs) mediante el pipeline DADA2 (v1.6.0; 
Callahan et al., 2016) que adicionalmente detecta y elimina secuencias quimera; 
corrigió el error por etiquetado de muestra (tag-switching) mediante la detección de 
las secuencias de la comunidad SynMock y el cálculo del indice de contaminación 
lateral y realizó un curado de ASVs erróneas mediante algoritmo LULU (v0.1.0; 
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Fr¢slev et al., 2017). Finalmente, el pipeline aplicado asignó taxonomía a las ASVs 
constituidas contrastando sus secuencias con la base de datos UNITE mediante un 
abordaje híbrido que combina alineamiento global y los algoritmos UTAX (RC Ed- 
gar, http://drive5.com/usearch/manual9.2/cmd_utax.html) y SINTAX (Edgar, 2016). 
La información de salida del pipeline aplicado se curó manualmente siguiendo las 
recomendaciones de Brown et al. (2015): las ASVs representadas por menos de 10 
lecturas, fueron eliminadas del conjunto de datos y se excluyeron de los análisis. 

Se extrajeron del set de datos las ASVs asignadas a los órdenes Polyporales e 
Hymenochaetales y su taxonomía se revisó manualmente. Cada secuencia se consul- 
tó en la base de datos GenBank (https://blast.ncbi.nlm.nih.gov/ a la que se accedió 
en noviembre de 2021) utilizando el algoritmo BLASTn (Altschul et al., 1997). La 
identidad del 90 al 95 % se consideró como una coincidencia a nivel de género, la 
identidad por encima del 97 % se consideró una coincidencia a nivel de especie. 
Cuando hubo conflicto entre las taxonomías de las múltiples secuencias consultadas, 
se asignó el último ancestro común. 


Análisis de datos 


Las diferencias de riqueza por muestra entre tratamientos se testearon mediante 
prueba no paramétrica de Kruskal-Wallis. Las diferencias de riqueza total entre tra- 
tamientos se testearon mediante test de Chi-Cuadrado de Pearson (1900). Diferencias 
de estructura de comunidad se testearon mediante prMANOVA (Anderson, 2017) 
utilizando un modelo que incorporara todas las variables de interés aplicado sobre 
matriz de distancias construidas con el índice de Raup-Crick pues es más adecuado 
para operar con datos binarios en sets de datos con muchos ceros (Raup y Crick, 
1979). Los análisis estadísticos y gráficos fueron realizados con el software RStudio 
(v1.1.456, http://www.rstudio.com/) con los paquetes biomformat (v1.8.0; McMurdie 
y Paulson, 2016), phyloseq (v1.24.2; McMurdie y Holmes, 2013), ggplot2 (v3.3.3; 
Wickham, 2016), vegan (v2.5.7; Oksanen et al., 2020) y betapart (v1.5.6; Baselga et 
al., 2022). 


RESULTADOS 


El experimento de secuenciación resultó en 38 millones de lecturas crudas y una pro- 
fundidad media de 133.855 lecturas por muestra. La curva de rarefacción se aproximó 
a una asíntota cuando se consideraron los datos seleccionados manualmente (datos 
no incluidos), lo que indica un muestreo satisfactorio. 

El pipeline aplicado arrojó 5726 ASVs informados, de los cuales se identificaron 
359 (6 %) a nivel de especie, 365 (11 %) a nivel de género, 280 (5 %) a nivel de familia, 
447 (8 %) a nivel de orden, 318 (6 %) a nivel de clase, 1663 (29%) a nivel de phylum 
y 2024 (35%) a nivel de reino. Se obtuvieron un total de 36 ASVs pertenecientes a 
los órdenes Polyporales e Hymenochaetales, detectados en 51 muestras de madera 
de Nothofagus. La lista de estos ASVs con sus taxonomías curadas y sus frecuencias 
de ocurrencia para cada variable evaluada se muestra en la Tabla 1. La riqueza por 
muestra fue mayor en N. dombeyi que en N. pumilio (Fig. 2). 
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Tabla 1. Taxones degradadores de los órdenes Polyporales e Hymenochaetales detectados en albura 
de individuos vivos de Nothofagus y sus frecuencias de ocurrencia para cada nivel de las variables 
evaluadas. 


Table 1. Degrading taxa of the orders Polyporales and Hymenochaetales detected in the sapwood 
of living individuals of Nothofagus and their frequencies of occurrence for each level of the evalu- 
ated variables. 


N. dombeyi N. pumilio 
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Nd: N. dombeyi. Np: N. pumilio. St: sintomático. At: asintomático. F: fuste. R: raíz. O: otoño. P: primavera. AC: sitio 
Cerro Alto el Petiso rodal de N. dombeyi. K: sitio Lago Kruger. ME: sitio Lago Menéndez. MS: sitio Rio Alerce. AL: sitio 
Cerro Alto el Petiso rodal de N. pumilio. D: sitio Cerro El Dedal. Ri: sitio Cerro El Riscoso. 
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Fig. 2. Diferencias en la riqueza por muestra de taxones de los órdenes Polyporales e Hymenochae- 
tales detectados en muestras de albura de N. dombeyi (Nd) y N. pumilio (Np) (prueba de Kruskal- 
Wallis). 


Fig. 2. Differences in the richness per sample of taxa of the orders Polyporales and Hymenochaetales 
detected in sapwood samples of N. dombeyi (Nd) and N. pumilio (Np) (Kruskal-Wallis test). 
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Independientemente del hospedante, el estado sanitario y el sitio, especies del 
género Postia fueron las más detectadas, seguidas por especies de Ganoderma y Phlebia 
(Fig. 3 y 4). El mismo patrón fue observado cuando las abundancias absolutas de 
lecturas fueron consideradas. 

Entre todos los sitios muestreados, hubo dos rodales de N. dombeyi que exhibie- 
ron valores más altos de riqueza total y riqueza por muestra: Río Alerce (13 ASVs) 
y Lago Kruger (15 ASVs). El sitio con valores más bajos de riqueza fúngica de los 
órdenes Polyporales e Hymenochaetales fue el rodal de N. pumilio en el Cerro El 
Dedal (1 ASV) (Tabla 2). 

La riqueza total de ASVs resultó diferente entre las especies de Nothofagus, 
siendo N. dombeyi la especie más rica y que presentó, ademas, mayores frecuencias de 
ocurrencia de taxones de los órdenes Polyporales e Hymenochaetales. N. pumilio ex- 
hibe una mayor proporción de presencias únicas que N. dombeyi (Tabla 2, Fig. 4). 

El análisis multivariado de la comunidad fúngica para todo el conjunto de datos, 
permite inferir que las mayores diferencias a nivel de la estructura de la comunidad 
se deben al compartimento vegetal (PERMANOVA: 3,27% de la varianza) y, el es- 
tado de salud (PERMANOVA: 2,96% de la varianza) (Tabla 2, Fig. 5). Los árboles 
de N. dombeyi y de N. pumilio decaídos mostraron una menor frecuencia fúngica de 
los órdenes Polyporales e Hymenochaetales que los sanos (Fig. 4) resultando en una 
diferencia de estructura significativa entre estados sanitarios en N. dombeyi (Tabla 2, 
Fig. 5). Las comunidades de N. dombeyi deterioradas se caracterizan por un conjunto 
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Fig. 3. Frecuencia de cada taxón de los órdenes Polyporales e Hymenochaetales detectados en 
muestras de albura de Nothofagus vivos. 


Fig. 3. Frequency of each taxon of the Polyporales and Hymenochaetales orders detected in sapwood 
samples of living Nothofagus. 
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Fig. 4. Frecuencias absolutas de taxones de los órdenes Polyporales e Hymenochaetales determinados 
usando una similitud de secuencia del 97 % con secuencias de referencia en la base de datos UNITE. 
Comparaciones para cada nivel de las variables estudiadas (hospedador, estado sanitario, órgano, 
estación y sitio) para todo el set de datos y para cada especie de Nothofagus por separado. Nd: N. 
dombeyi. Np: N. pumilio. St: sintomático. At: asintomático. F: fuste. R: raíz. O: otoño. P: primavera. 
AC: sitio N. dombeyi Cerro Alto el Petiso. K: Lago Kruger. ME: Lago Menéndez. MS: Río Alerce. AL: 
sitio N. pumilio Cerro Alto el Petiso. D: Cerro El Dedal. Ri: Cerro El Riscoso. 


Fig. 4. Absolute frequencies of taxa in the orders Polyporales and Hymenochaetales determined 
using 97% sequence similarity to reference sequences in the UNITE database. Comparisons for each 
level of the variables studied (host, health status, organ, season and site) for the entire data set 
and for each Nothofagus species separately. Nd: N. dombeyi. Np: N. pumilio. St: symptomatic. At: 
asymptomatic. F: stem. A: root. Or: autumn. Q: spring. AC: Alto el Petiso Mt. N. dombeyi site. K: 
Lake Kruger. ME: Lake Menéndez. MS: Alerce River. AL: Alto el Petiso Mt. N. pumilio site. D: El Dedal 
Mt. Ri: El Riscoso Mt. 


de 13 taxones únicamente detectados en individuos sintomáticos, mientras que las 
comunidades de árboles sanos podrían caracterizarse por un conjunto de 10 taxones 
únicamente detectados en individuos visiblemente saludables (Tabla 2, Fig. 6). 

La riqueza total de hongos degradadores en raíz fue mayor que en fuste en N. 
pumilio (Tabla 2). Las comunidades de fuste y raíz de N. dombeyi se diferenciaron en 
estructura presentando riquezas similares (Tabla 2, Fig. 5). 
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Tabla 2. Diversidad de taxones de los órdenes Polyporales e Hymanochaetales en albura de Notho- 
fagus vivos. Riqueza de cada nivel de las variables estudiadas, diversidad beta entre ellos y efecto 
de las variables estudiadas sobre la estructura comunitaria. 


Table 2. Diversity of taxa of the orders Polyporales and Hymanochaetales in living Nothofagus sap- 
wood. Richness of each level of the variables studied, beta diversity between them and effect of the 
variables studied on the community structure. 


Diversidad Alfa Diversidad Beta 
Chi-cuadrado PERMANOVA 
riqueza p-valor p-valor R2 Sorensen turnover  nestedness 
hospedador Na 30 ie 0,0184 0,707 0,500 0,207 
Np 123 
estado sanitario a a4 0,0296 0,707 0,700 0,007 
At 21 él 
órgano r 13 0,0327 
R 22 ig 
estación E a2 0,0193 
P 24 
estado sanitario a 13 0,0497 0,758 0,733 0,024 
At 16 + 
á F 18 
órgano 0,0566 0,758 0,750 0,008 
R 17 ia 
A ee 0) 20 
N. dombeyi estación P 15 0,0285 0,657 0,625 0,032 
AC 7 
sitio K a3 0,0817 
ME 7 
MS 13 
ates ts St 6 
estado sanitario al 7 0,0833 0,846 0,833 0,130 
F F 3 
órgano i 0,0833 0,846 0,667 0,179 
R 10 
N. pumilio estacion > x 0,0833 0,846 0,750 0,096 
AL 7 
sitio D 1 0,1667 
Ri 2 


Nd: N. dombeyi. Np: N. pumilio. St: sintomático. At: asintomático. F: fuste. R: raíz. O: otoño. P: primavera. AC: sitio 
Cerro Alto el Petiso rodal de N. dombeyi. K: sitio Lago Kruger. ME: sitio Lago Menéndez. MS: sitio Río Alerce. AL: sitio 
Cerro Alto el Petiso rodal de N. pumilio. D: sitio Cerro El Dedal. Ri: sitio Cerro El Riscoso. 


DISCUSIÓN 


Los estudios con abordajes metabarcoding proliferaron las últimas dos décadas y evi- 
denciaron la existencia de una gran diversidad oculta, aún desconocida, que podía 
involucrar inclusive, linajes enteros (Grossart et al., 2016; Hibbett et al., 2016). En 
los bosques templados del hemisferio sur y particularmente para Patagonia, todavía 
son muy escasos los estudios con estas nuevas tecnologías en hongos, además de 
que se sabe muy poco sobre la interacción y las funciones de dicha micobiota en 
el decaimiento de los bosques, la descomposición de la madera y la muerte de los 
árboles (Baldrian, 2016). 

El presente trabajo amplía los hallazgos de estudios anteriores y aumenta nuestra 
comprensión de la biología de los hongos que causan la descomposición de la madera 
en los bosques nativos de Patagonia, Argentina. De manera más general, contribuye 
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Fig. 5. Escalado dimensional no métrico (NMDS) reflejando las diferencias de estructura de las co- 
munidades entre órganos (izq.) y estados sanitarios (der.) para todo el set de datos (arriba) y para 


N. dombeyi (abajo). 
Fig. 5. Non-metric dimensional scaling (NMDS) showing differences in community structure between 
organs (left) and health status (right) for the entire data set (top) and for N. dombeyi (bottom). 
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Fig. 6. Taxones de los órdenes Polyporales e Hymenochaetales exclusivos (en color) y compartidos 
(negro) de cada especie de Nothofagus y de individuos sintomáticos (rojo) y asintomáticos (azul). 
Fig. 6. Exclusive (in color) and shared (black) taxa of the Polyporales and Hymenochaetales orders 
of each Nothofagus species and of symptomatic (red) and asymptomatic (blue) individuals. 
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al conocimiento más amplio de este grupo a nivel internacional, particularmente 
en el hemisferio sur, donde se encuentran varias especies de basidiomicetos que no 
están presentes al norte del ecuador. En este estudio, describimos las comunidades 
de hongos causantes de pudriciones de los órdenes Polyporales e Hymenochaetales 
de las dos principales especies de Nothofagus de los bosques patagónicos andinos del 
norte, complementando los hallazgos basados en cultivo con resultados de enfoques 
libres de cultivo (Molina et al., 2020). Se amplificó la región ITS1 y se envió a NGS, 
generando secuencias aleatoriamente. Las longitudes de lectura obtenidas fueron 
suficientes para discriminar unidades taxonómicas moleculares operativas en el nivel 
de género al nivel de similitud del 97% dentro de estos grupos fúngicos. 

Hemos identificado un mayor número de géneros de hongos degradadores en N. 
dombeyi que en N. pumilio. El efecto del hospedante en las comunidades de hongos 
ha sido informado como fuerte, bajo e incluso no significativo en diferentes estudios 
(Bahram et al., 2022; Kiingas et al., 2020; Rim et al., 2021; U’Ren et al., 2019, respec- 
tivamente). La divergencia de las comunidades fúngicas de N. dombeyi y N. pumilio 
se ha discutido como un reflejo de las distribuciones altitudinales no adyacentes en 
el norte de la Patagonia, la disponibilidad de agua y la edad del rodal (Molina et al., 
2020). Mientras que en Nueva Zelanda vincularon la jerarquía de similitud con las 
distancias filogenéticas entre las tres especies de Nothofagus evaluadas (Johnston et 
al., 2012, 2017). 

La composición, riqueza y diversidad de los hongos degradadores se rigen por 
muchos factores, tanto externos como intrínsecos. Algunos basidiomicetos son inhe- 
rentemente más vigorosos y capaces de superar a otros, y las especies pueden emplear 
naturalmente diferentes estrategias de ocupación. La capacidad de colonización de 
diferentes hospedantes podría estar relacionada con propiedades particulares del 
duramen o la albura. Los hongos de la pudrición de la madera difieren de forma 
innata en su capacidad para producir enzimas capaces de desnaturalizar la celulosa 
y la lignina (Hood et al., 2019). 

Especies de Postia y Ganoderma fueron los taxones inequívocamente dominantes 
principalmente en árboles asintomáticos de N. dombeyi. Postia, causante de pudrición 
castaña, es uno de los géneros más grandes de Polyporales (Basidiomycota). Postia 
pelliculosa es una de las principales especies que pudren la madera de varias espe- 
cies de Nothofagus en los bosques andinos patagónicos y está muy extendida en el 
hemisferio sur, siendo registrados también en el este de África (Ryvarden y Johan- 
sen, 1980), Camerún (Douanla-Meli, 2007), Australia (Cunningham, 1965) y Nueva 
Zelanda (Buchanan y Hood, 1992). Ganoderma australe (Polyporales) importante en 
ambas especies de Nothofagus, pero acompañada por un amplio rango de otras espe- 
cies de Basidiomycetes, ha sido reportada como causante de pudrición blanca en la 
región (Rajchenberg, 2006). Las especies de Ganoderma han sido registradas como 
degradadores y patógenos en otras especies de Nothofagus en Nueva Zelanda y en 
muchas especies vegetales alrededor del globo (Hood et al., 2019). 

La estación del año fue otro factor que afectó la estructura de la comunidad. El 
muestreo de primavera representa el comienzo de esta temporada y refleja el efecto 
de los períodos más fríos, que también son los más húmedos. Por el contrario, al final 
de la temporada de crecimiento (muestreo de otoño), podemos observar el efecto del 


168 L. Molina, M. B. Pildain: Hongos degradadores de la madera en bosques Andino-Patagónicos 


período más cálido y seco siguiendo un patrón retrasado como en los suelos (Kunkel 
et al., 2016; VoYiskova et al., 2014). En primavera encontramos una mayor riqueza de 
taxones y una mayor abundancia absoluta que en otoño y esto ocurre principalmente 
en N. pumilio. Estos patrones, acompasados a la estacionalidad, ya se han observado 
para la micobiota asociada a raíces (Wang et al., 2021). De esta manera, la comunidad 
de hongos degradadores de N. pumilio presentaría una mayor susceptibilidad a la 
estacionalidad y a los cambios de temperatura, en comparación con el micobioma 
de N. dombeyi, que se distribuye en ambientes menos extremos. 

Por otro lado, la riqueza total de hongos degradadores en raíz fue mayor que en 
fuste en N. pumilio. Se conoce que la degradación en el suelo o en contacto con el 
suelo típicamente es más rápida, reflejando la presencia de mejores condiciones de 
humedad acompañada con una abundancia de organismos con diferentes capacida- 
des de degradación de la madera. En contraste, el ambiente sobre el suelo donde se 
encuentra el fuste, es más desafiante, con menos organismos presentes y más seco, 
lo cual puede enlentecer el proceso de degradación y provocar cambios sustanciales 
en la estructura de las comunidades presentes (Goodell, 2020). 

Reflejos de los patrones de diversidad y estructura de hongos degradadores uti- 
lizando nuevas herramientas son estudios emergentes, especialmente para bosques 
templados del hemisferio sur. Las especies de árboles analizadas necesitan de aportes 
en relación a las funciones, condiciones y procesos tomado como un holobioma en 
que los hongos son como ejes centrales de cambios de comunidades que repercuten 
en el funcionamiento del organismo. Es necesario ampliar el conocimiento sobre 
hongos claves en la sucesión para comprender los motores de diversidad y estructura 
y la colonización de tejidos de estas especies forestales. 
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